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(i)  SUMMARY  OF  RESEARCH  COMPLETED  TO  DATE 
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1  .  A  survey  of  literature  concerning  ultrasonic  nondestructive  testing 
of  composite  materials.  This  in  included  in  Appendix  I. 

2.  A  theoretical  study  of  wave  phenomena  unique  to  the  (anisotropic) 
composite  materials  has  been  completed.  This  is  included  in  Appendix 
II. 

3.  A  theoretical  study  of  the  experimental  technique  under  development 
has  been  complteed  (Section  2) 

4.  A  very  accurate  computer-controlled  experimental  technique  to 

measure  wavespeed  in  a  composite  material  has  been  developed.  For 

* 

conventional  isotropic  materials  accuracies  of  two  parts  per  ten 
thousand  were  achieved;  for  highly  damaged  Graphite/epoxy  composite 
materials  accuracies  of  0.2$  have  been  achieved.  |C^^4-oc:  / 
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I.  INTRODUCTION 


It  is  well  known  that  components  made  of  composite  materials 
suffer  complex  damage  when  they  are  subjected  to  either  monotonic  or 
fatigue  loading.  We  submitted  a  three  year  proposal  to  develop 
ultrasonic  nondestructive  techniques  to  measure  the  damage  states  in 
fiber-reinforced  composite  materials  (Texas  A&M  Research  Foundation 
Proposal  No.  RF-8M-3^)  The  objective  of  this  Annual  Technical  Report  is 
to  summarize  the  progress  made  during  the  first  year  of  the  grant  period 
(AFOSR-8H-0066) . 

When  the  damage  occurs  it  has  two  effects  upon  the  propogation  of 
a  mechanical  wave  through  the  composite:  (1)  It  affects  the  stiffness, 
and  therefore,  the  speed  of  wave  propagation;  (2)  It  increases  the 
attenuation  of  the  wave.  Thus  by  measuring  the  speed  and  attenuation  as 
functions  of  frequency  at  various  known  levels  of  damage,  we  propose  to 
establish  the  necessary  calibration  curves  for  damage.  Conversely,  one 
can  measure  the  ultrasonic  parameters  and  estimate  the  extent  of  damage. 

In  Section  2  we  summarize  the  theoretical  results  obtained  to 
date.  In  section  3  we  describe  the  developilffent  of  the  experimental 
technique.  Section  4  is  concerned  with  the  implementation  of  the 
experimental  techniques  i.e.  with  the  measurement  of  phase  velocity  and 
attenuation  as  the  damage  progresses  in  graphite/epoxy  composite 
specimens 

2.  THEORETICAL  RESULTS 

In  this  section  we  describe  the  theoretical  foundations  of  the 
experiment. 

Consider  a  plate  of  a  linear  viscoelastic  material  sandwiched 
between  two  half-spaces  of  perfectly  elastic  materials.  Consider  a 

2 


finite-duration  pulse  incident  at  the  viscoelastic  plate.  Because  of 
the  mismatch  in  the  acoustic  impedance  of  the  three  materials  there  will 
result  an  infinite  series  of  reflected  and  transmitted  pulses.  The 
reflected  and  transmitted  fields  contain  information  about  the  acoustic 
properties  of  the  viscoelastic  material.  An  objective  of  this  section 
is  to  derive  expressions  for  the  phase  velocity  and  the  attenuation  of 
waves,  given  the  incident  and  the  reflected  (or  the  transmitted)  fields. 

In  view  of  the  complexity  of  the  total  problem  we  have  broken  it 
down  into  the  following  simpler  sub-problems.  We  will  first  assume  that 
the  plate  is  perfectly  elastic.  An  incident  pulse  of  arbitrary  shape 
will  be  considered.  The  most  general  expressions  for  the  reflected  and 
the  transmitted  fields  will  be  derived.  Next,  instead  of  a  pulse  of 
arbitrary  shape  we  will  consider  a  monochromatic  wavetrain.  For  a 
perfectly  elastic  material  the  wavenumber  k  is  real.  At  this  point  the 
viscoelasticity  of  the  composite  will  be  accounted  for  by  letting  the 
wavenumber  k  take  on  complex  values  (k=k^+ik2).  A  pulse  of  arbitrary 
shape  can  now  be  studied  in  terms  of  its  Fourier  components.  Finally 
from  this  general  analysis  we  will  derive  expressions  for  a  few 
particular  situations  which  will  be  used  in  the  laboratory. 

First,  consider  the  reflection  and  transmission  of  a  right-going 
displacement  wave  at  the  interface  of  two  elastic  media,  see  Fig.  1.  The 
equation  of  motion  is 

=  _1_  9^u  (1) 

2  2  2 
9x  C  at 

where  u  is  the  displacement  in  the  x-direction,  t  is  lim-  ind  C  is  the 

2 

speed  of  the  longitudinal  wave:  C  “(A+2u)/pi  wh'.re  K  and  u  are  Lame's 
constants  and  p  is  the  density.  Let  the  incident  pulse  be  given  by 


=  f(u)t-kx) 


I 


f(3)  s  0  for  s<0  (2) 

where  at  present  u  are  k  are  any  two  constants  connected  by  C=(ij/k. 
Later,  oj  and  k  will  be  identified  with  the  circular  frequency  and 
wavenumber  of  a  monochromatic  harmonic  wave.  The  boundary  conditions  at 
x=a  are  the  continuity  of  displacement  and  stress.  The  reflection  is 

given  by  (with  s=u)t+k  x) 

r  f 

u  =  R,_f(s-2k  a)  for  s>2k  a  (3) 

.  1 2  o  o 

0  for  3<2k  a 

^  o 

where  the  displacement  reflection  coefficient  is 

8,,- 

rcT^c  (“> 

o  o 

The  transmitted  field  is  given  by  (with  s=wt-kx) 
u'^-Jt^ s>(k^-k)a 

L  0  for  s<(k^*-k)a  (5) 

where  T^  ^  ^  ^ 

Next  consider  a  left-going  displacement  wave;  Fig.  1b 
u^"°=f(u)t+kx)  (7) 

f(s)=0  for  s<so 

where  is  some  arbitrary  constant. With  s=u)t-kx,  the  reflected  wave  is 

P  f 

u  =  R_,f(s+2ka)  for  s>s  -2ka 

0  for  s<s^-2ka  (8) 

With  s=a)t+k  x  the  transitted  wave  is 
o 

,  T„,  f  (s- (k  -k)a)  for  s>s  +(k  -k)a 
^21  o  o  o 

0  for  3<s  +(k  -k)a  (9) 

o  o 


^21  '^12  ”^21  ^  ”^1 


>  V .  .  .  .  .-a.  r>>^  -Ji  .-li  .'Ji 


We  next  consider  an  elastic  plate  sandwiched  between  two  elastic 


half-spaces.  A  convenient  way  of  following  the  various  reflections  is 
to  draw  a  t-x  diagram  of  a  point  of  constant  phase  in  the  incident 
pulse;  this  is  shown  in  Fig.  2.  This  method  is  also  known  as  ray 
tracing.  Let  the  displacement  field  along  the  incident  ray  (ray  1  )  be 
given  by 


uinc_f(^t_[^  x) 
o 


(11  ) 


f(s)=0  for  s<o 

The  reflection  (ray  2)  and  transmission  (ray  3)  at  x=a  are  given  by  (in 
the  following  the  superscript  refer  to  the  corresponding  ray  number): 


u  = 


0  for  s<s. 


(12) 


where  s=(ijt-k  x  ahd  s„=2k  a,  and 
o  2  o 


s>S2 


for  3<s, 


(13) 


where  s=a)t-kx,  s.  =  (k  -k)a 
3  o 


and 


„  p  C-oC  .  „ 

Ri2  =  VVt  12 

0  C  +pC 
o  o 


2p  C 
o  o 

p  C  +pC 
o  o 


(14) 


Ray  3  now  suffers  a  reflection  (ray  5)  and  a  transmission  (ray  4)  at 


x=b. 


Ti2fT23f(s-Si^)  for  s>s^ 

0  for  s<Sj^, 


with  s=u)t-k-x  and  s.=k  a+k(b-a)-k.b 
3  4  o  3 


T,  _R.,_f  (s-s^)  for  s>s^ 
12^3  b  b 

0  for  s<s^ 


with  .s=ait+kx  and  s^=k  a+k(?b-a)  and 
5  o 


(15) 


(16) 
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R  =  ^3^3  and  T  =  2pC 

oC+p^C^  oC+P^C^ 


In  a 
a  Ion 
foil 


similar  manner  one  can  write  down  the  expressions 
^  various  rays.  The  total  reflected  field  is 


owing  rays  (with  h=b-a  and  s=ujt+k^x); 


kh 


u"  =  R,  .,f  (s-3,  ) ;  S,,=--Zk  a 
12  2  2  O 

^  2^23'^21 

u  =’^1  2''^23^21  ^23^21 


for  displacement 
the  sum  of  the 


The  sum  may  be  written  as 

”  m-1 


where  s=ajt+k  x  and  s  =2k  a+m2kh. 
o  mo 

In  (17)  m  is  the  number  of  complete  round  trips  taken  by  the  wave  across 
the  plate  thickness  h 

In  an  exactly  analogous  manner  one  can  write  the  expression  for 
the  total  transmitted  field.  With  s*=a)t-k^x 

u^=T^2'^23^^"'^^’ 

u^=T^  ^23^ ’  Sg=a(k^-k2)+h(3k-k2) 

1 2  2 

u  =T.|  2^f^23^21  ^  T^^f  (s-s^  2^ !  o"3(k^-k2)+h(5k-k2) 


The  total  transmitted  field  rr.av  be  written  as 


00 


“  ■’■l2’'23  ^  '’’23"21>' 

m=(J 

s  =a(k  -k.  ) +:i[  (  2n'i  +  1  )k-k.  ' 
ni  o  5  3 


1  3; 


It  is  emphasised  that  in  the  foregoint^  tiie  plate  is  assurr:oc  to 
behave  in  a  penfeetly  elastic  manner'.  It  is  noted  that  with  a  sligiit 
ciiange  of  detail  the  precoeding  analysis  is  equally  valid  for  an 
incident  shear  wave. 

Water  Immersion  Method 

We  now  consider  the  mechanical  response  of  a  linear  viscoelastic 
plate.  If  the  attennation  is  small  this  can  be  most  easily 

accomplished  by  letting  k  in  the  expression  for  a  time-harmonic  wave, 
kx)^  take  on  complex  values.  Consider  a  basic  experimental  set-up 
shown  in  Fig  3(b).  The  piezoelectric  transducer  (T)  and  the  composite 
specimen  are  acoustically  compled  through  water.  An  electrical  signal 
is  applied  to  the  transducer  which  launches  a  longitudinal  wave  in  the 
direction  of  the  composite.  The  first  order  of  business  is  to 
characterize  this  pulse.  To  this  end  we  measure  the  reflection  from  a 
homogeneous  medium  of  known  properties  (p  and  c);  Fig  3(a)  i  Let  the 


incident  field  be 
inc 


u 


f  (oit-k  x). 
o  o 


(19) 


The  reflection  is  given  by 

p  P  C  "*pC 

u  =Rf  (ut+k  x-2k  a):  R=  o  o  (20) 

°  °  °  PqCo+pC 

Let  the  reflected  signal  as  sensed  by  the  transducer  be  labelled  f(t), 
then 

f  (t)=f^(ait-2k^a)  (21) 

Define 


F* 


"L 


(t)  dt 


and  F  *(0))  -vni  f  (<*3t)  dt. 


(22) 

(23) 
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(24) 


! 


j- 


V 

sJ 


V 


.1 


then  F  »(a))=  F»(u)) 

o 


Note  that  F*(w)  is  the  Fourier  transform  of  a  measured  quantity,  f(t), 
from  which,  via  eq.  (24),  we  can  deduce  the  Fourier  transform  of  the 
signal  launched  by  the  transducer.  Now  consider  a  wave  travelling  in 


the  positive  x-direction. 

inc,  .  ,  _L_  r„  * 

u  ( X , t )  V^tTj 


,  .  i((ut-k  x) 
(tjj)  e  o  dcj 


(25) 


Here  k  is  pure  real.  In  comparison  with  eq  (19) 


1  r 

f  (tot-k  x)='^TIf 
o  0  ^  J  o 

:e 

1 


+  <» 


i(wt-k  x)  . 
*(u)  e  o  dw 


(26) 


Hence 

r(3)WJ7/p„  *(0))  e^^  du  (27) 

—CO 

The  reflected  field  is  given  by  eq.  (17).  Setting  •x=0  (since  the 
transducer  is  located  at  x=o)  the  total  reflected  field  is  given  by 
u^(o,t)^\5T|f’o*(w)da) 


-i2k  a  ^ 
e  o  + 


^6  g-i{2kQa+2mkh}J 


(28) 


m=l' 


Here  ^^^21  ”^21  ^21 


2Cm-1) 


Let  u^(o,t)  be  labelled  as  g(t)  and  let 


I  “ 

G»(a))=V?^g(t)  e  dt. 


..  r„  -i2mkh  2ik  a  „ 

then  2,6  e  e 

^  m  F*  12 


(29) 


By  observing  that 


1+Z+Z+  + . =  JL  for  1Z1<1 

1-Z 

the  foregoing  expression  can  be  recast  in  the  form 


„  2  -i2kh 
R21  e  =6 


H-B 


T  T 
12  21 


[«  1:  -  «12] 


(30) 
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I, 


This  is  the  key  result  of  this  section.  By  measuring  F*  and  G*  one  can 
extract  the  complex- valued  k(!i))=k^ +ik2  from  eqs.  (30) 

Measurement  of  k(u))  from  the  transmitted  field 

An  alternative  method  of  measuring  k(m)  is  from  the  transmitted 
field.  With  reference  to  Fig.  3  a  matched  ultrasonic  transducer  is  used 
as  a  receiver  (R).  The  reference  signal  f(t),  is  obtained  by  removing 
the  specimen  from  the  wave  path.  The  incident  field  is  given  by  f^  (tot- 
kx).  Therefore 

f(t)=fo(a)t-kl)  (31) 

Let  F*(u))  and  F^*(!jj)  be  defined  as  in  equations  (22)  and  (23), 
respectively,  then 

F*(io)=e^^°^F»(oj)  (32) 

Let  g(t)  be  the  total  transmitted  field,  ray  4,8,12 . “>.  Let 

1  +” 

G*(a))=V2i“  /g(t)e"^‘^^  dt  (33) 

Then  from  eq.  (18) 

G»(a))=T  T..  F»((u)  e"^‘^‘^.  I  (34) 

21  ^=0 

where  z=[E^.  e 

observing  that  ^  eq  (34)  reduces  to 

m=0  1-Z 


-ih(k-k  ) 
o 

e _ 

1  R  2  -i2kh 
^■‘^21  e 


G*(u) 


T12T21FMU.) 


F*(a))  and  G*((d)  are  measured  quantities,  therefore,  eq.  (35)  is  viewed 
as  an  equation  to  evaluate  k((jo). 

For  very  thin  composites  it  is  frequently  impossible  to  separate 
the  various  rays  (Fig.  2)  in  the  time  domain.  The  preceding  discussion 
deals  with  that  general  case.  Occasionally  it  is  possible  to  isolate 


various  rays.  when  tnat  happens  certain  simplifications  result.  In 

particular  let  us  suppose  that  it  is  possible  to  isolate  the  front 

surface  reflection  (ray  .-1)  from  the  remaining  reflected  field.  Then  it 

is  more  convenient  to  use  ray  2  as  the  reference  signal  i.e. 

f  (t)=R^  (36) 

and  F  »(w)=_l_e^"'^o'^  F*(u))  (37) 

°  ^12 

Define  g(t)  as  the  total  reflected  field  after  ray  2  i.e.  ray  6, 


Then  the  results  of  a  calculation  similar  to  the  one 


described  in  the  foregoing  yield: 


i2kh  ^  , 
e  =  H-  ^  , 


where  6=  ^21^12  G^-  , 

T  T  F* 

12  21 

which  is  the  equation  for  calculating  kCu). 

Finally,  we  treat  the  case  when  the  damage  is  large  and  therefore, 
the  damping  is  high.  When  that  happens  Rays  10,  1^, ...  .become  of 

insignificant  amplitude.  The  wavespeed  and  attenuation  are  best 
measured  by  convoling  Ray  2  with  Ray  6.  If  f(t)  is  the  transducer 
response  to  Ray  2,  g(t)  the  response  to  Ray  6,  F*(u))  and  G*(a))  are  the 
Inverse  Fourier  Transforms  of  f(t)  and  g(t),  respectively,  it  can  be 


readily  shown  that 

G*M  =  T,„T  e 
F*(u)) 


-i2kh 


Direct  Contact  Method  In  addition  to  the  water  immersion  method  we  have 
used  a  direct  contact  technique  in  which  the  transducer  and  the  specimen 
are  acoustically  coupled  with  a  block  of  an  isotropic  homogeneous 
material  (e.g.  aluminum)  rather  than  water.  A  schematic  arrangement  is 


shown  in  Fig.  4.  Within  this  arrangement  we  will  discuss  three  cases 


which  were  found  useful  in  the  laboratory. 


Case  1 . 


Reference  Signal:  Reflection  from  the  back  surface  of  the  delay  rod 

with  the  specimen  removed. 

Reflected  field:  Rays  2,6,10,14 . “ 

u  =f  (ut-k  x) 
o  o 

f(t)=f  (a)t-2k  a)  (40) 

0  o 

F*(u))=e“^^'^  ^Fo*(w)  (41) 

oo  m-1 

I  ^2) 

m=l 

where  s  =2k  a+2mkh 

mo  oo  m—  1 

m=l 

After  some  algebraic  manipulation  we  get 


R2ie"^^‘^'^=8/(1+8) 

1 

where  6».— iJ: _ 

^12^21 


(G*  -  R.„) 

F* 


which  are  the  desired  equations  for  calculating  the  complex-valued  k(a)) 
Case  2 

Reference  Signal:  Ray  2 

Reflected  Field:  Rays  6,10,14,....® 

u^’^°=fo(u)t-kx) 

f(t)=R, -fo(a)t-2k  a)  (46] 

1 2  o 

F*(a))=R^2  e"^^‘^o^Fo»(a))  (47] 

oo  m-1 

g(t)=Ti2T2i  I  R2^fo(a)t-s  )  (48] 

where  s  =2k  a+2mkh 
m  o 

oo 

r*(  ^  TioToi  ^  „  -i { 2koa+2mkh }  \  ciin' 

G*(a))=__ti_£l  I  Rp^  e  Fo*((jj)  (49. 

m=l‘^ 


After  some  algebraic  manipulation,  it  can  be  shown  that 

giSkh^j^^^  (^F*  +  g»)/G»; 
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*»*? 


where  a  2  "^21^^12^21 
Case  3 

Reference  Signal:  Ray  6 
Reflected  field:  Rays  6,10,1^1 . 


When  the  excitation  voltage  is  applied  to  the  transmitting  transducer, 
it  saturates  the  preamplifier  in  certain  situations.  Frequently,  there 
is  enough  time  for  the  preamplifier  to  have  a  complete  recovery  before 
the  first  reflection,  ray  2,  arrives;  sometimes  there  is  not  enough 

time.  This  calculation  is  done  for  these  situations. 

inc  «  /  .  ,  \  4 

u  =f-(a)t-k  x) 
o  o 


f  (t)=T,  _T-,fo(a)t-s.) ;  s^=2kh+2k  a 

12  21  D  0  O 

(51  ) 

F»(a))=T^2'^21 

(52) 

»  m-1 

g(t)«Ti2T2i  I  R2^fo^“''"®o)‘»  s^=2k^a+2mkh 

ID=1 

(53) 

«  m-1 

G»(u))=T^2'^21  ^  ^21  Fo»(a)) 

ra=l 

(54) 

R21  e"^^‘'‘"-(GtF*)/F* 

(55) 
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3.  TECHNIQUE  DEVELOPMENT* 

Keeping  in  view  the  tremendous  speed  and  reliability  we  can  achieve 
by  using  computers  for  collection  and  analysis  of  data,  equations 
developed  in  section  2  were  interpreted  in  a  way  most  suitable  for 
computer  programming.  The  technique  was  developed  such  that  the 
meaurements  could  be  made  with  least  amount  of  operator  interference. 

The  technique  thus  developed  was  first  applied  to  aluminum  samples 
as  a  calibration  measure.  The  reason  being,  the  ultrasonic  waves  are 
weakly  dispersive  in  aluminum  and  the  wavespeed  is  known  apriori.  The 
transducer  was  kept  in  direct  contact  with  the  specimen,  coupled 
acoustically  with  a  thin  layer  of  mineral  oil.  The  ot*her  end  of  the 
specimen  was  free,  in  contact  with  air,  see  Fig.  4. 

When  the  transducer  is  excited  by  a  short  duration  pulse,  it 
launches  a  wave  in  the  specimen  and  this  wave  reflects  back  and  forth 
between  the  parallel  surfaces  of  the  specimen  giving  a  signal  shown  in 
Fig.  5.  The  signal  consists  of  pulses  and  for  later  discussion  the 
pulses  have  been  numbered.  It  is  observed  from  this  figure  that  there 
are  regions  between  pulses  where  the  signal  level  is  very  low.  We  can 
clearly  define  the  end  of  one  pulse  and  the  start  of  the  next  pulse. 
This  separation  of  pulses  depends  on  the  excitation  frequency,  sample 
thickness,  wavespeed  and  dispersion  in  the  specimen.  As  a  first  step  the 
measurement  technique  is  developed  on  thick  specimen  where  this 
separation  is  possible.  The  reason  why  we  are  emphasizing  on  pulse 
width  will  become  clear  during  subsequent  discussions.  In  the  case  of 
composites  we  will  not  be  able  to  avail  this  luxury  and  there  we  will 
use  the  equations  in  such  a  way  that  measurement  will  still  be  possible. 


*Prepared  by  Vinay  Dayal  (Ph.D.  Student) 


The  equation  derived  in  section  ?.  for  the  transducer  in  contact 
with  specimen  is  reproduced  here  for  the  continuity  of  discussion. 

G»(w)-FMw)  =  R  ,  (55) 

F»(w) 

for  aluminum  k  is  real.  We  now  define  t^  as  the  round  trip  time  (t 

o  o 

=  2h/C),  and  the  circular  frequency  =kc,  then 

G»(w)=F*(’j.')  [1 (56) 

here  F*(w)  is  the  Fourier  Transform  (FFT)  of  pulse  I  of  Fig.  5. 
(Ray  6  of  Fig.  2)  and  G*(u))  is  the  combined  FFT  of  I  and  II  (Ray  6  & 
10).  Both  F(a))  =  lF*(aj)|  and  G(u))  =  |G*(a))j  are  plotted  in  Fig.  6.  It  can 
be  seen  that  F((d)  is  an  inverted  bell  shaped  curve  and  is  essentially 
the  transducer  response.  While  G(io)  exhibits  harmonics  due  to 
constructive  and  destructive  interference  of  the  first  and  second  pulse, 
which  is  superimposed  on  F(u). 

If  Af  is  the  difference  in  frequency  between  the  two  successive 
peaks,  it  can  be  readily  shown  that  the  wavespeed  is 

C=2hAf  (57) 

From  Fig.  6,  average  peak  spacing  on  frequency  scale  is  measured 
and  wavespeed  is  calculated. 

Another  technique  of  measuring  c  is  to  recast  equation  (56)  in  the 

form 

P*(^)  -1=R^,  (58) 

F*(a0 

and  plot  1^,-1 1  vs  frequency.  The  results  are  shown  in  Fig.  7. 

F* 
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+  2A  COS(a)t 

F(u)) 
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Here  we  have  deconvolved  the  syste.r,  response  from  the  interference 
Pattern  of  two  pulses.  With  Af  now  interpreted  as  the  frequency 
spacing  between  two  successive  zero  crossings,  equation  (57)  still  holds 

for  cdlculat ing  C- 

These  two  methods  were  used  on  aluminum  specimen  and  wavespeed  was 
found  to  be  accurate  to  ±  0.5?.  The  main  drawback  of  these  techniques 
is  that  the  peak  spacing  had  to  be  measured  manually  and  also  in  cases 
where  the  number  of  peaks  becomes  small,  it  will  be  difficult  to  get 
accurate  peak  spacing.  As  stated  earlier  our  objective  in  this  research 
effort  is  to  develop  a  fully  automated  computer  controlled  measurement 
system.  We  proceed  ahead  and  recast  equation  (50)  in  the  form 

G*(aj)  G*(a))-F*(u0  E  S*(aj)  ^  R., 

F*(oj)  F*(a))  F*(w) 


where  S*(w)  is  the  Fourier  Transform  of  the  second  pulse  (Fig.  5). 
Both  S*(w)  and  F*(w)  are  complex  numbers, equation  (59)  can  be  written  as 

S*(a))  ^ 

=  |Mlei'^  (60) 

where  if  is  the  phase  difference  between  the  phase  of  second  pulse 

and  the  phase  of  the  first  pulse.  Since  we  are  interested  in  the 
difference  in  phase  we  do  not  care  about  the  absolute  phase  and,  taking 
the  first  point  on  the  signal  as  the  .•J)=0  point,  calculate  the  phase 
difference  between  the  two  pulses  at  each  digitizing  point. 

Comparing  equation  (59)  and  (60)  we  get 


Thus  6  =  't  =2v  f  2h 


HERTZ 


or  <{>  =  (A-h)  f 
c 

Now  if  we  plot  the  phase  difference  ({>  vs  the  frequency  we  should 
obtain  a  straight  line  with  a  slope  of  AttN/c.  Fig.  8.  is  the  phase  vs 
frequency  plot  for  an  aluminum  specimen.  It  is  seen  that  over  the  range 
of  about  4.5-10  MHz  a  fairly  straight  line  is  obtained.  This  portion  of 
the  curve  is  expanded  for  clarity  in  Fig.  9.  This  straight  line  section 
of  the  curve  is  related  to  the  frequency  response  of  the  transducer  and 
will  be  discussed  later.  A  least  square  fit  line  is  passed  through  these 
points  and  the  slope  is  computed  from  which  the  wavespeed  is  calculated. 
All  these  calculations  are  performed  on  a  computer  ‘with  no  human 
intervention. 

This  form  of  equation  (60)  has  been  used  in  the  development  df  the 
computer  algorithm  and  the  measurements.  Now  we  will  discuss  some 
factors  which  affect  the  accuracy  of  our  measurements. 

The  first  factor  is  the  digitizing  interval  for  the  signal.  The 
analog  signal  in  time  domain  is  digitized  by  the  oscilloscope  at  fixed 
intervals.  We  would  like  the  digitization  to  be  done  at  such  an 
interval  that  no  part  of  the  signal  is  lost.  If  the  digitizing  interval 
is  too  large  it  would  result  in  the  loss  of  high-frequency  components  of 
the  waveform.  This  effect  was  studied  by  digitizing  a  signal  at  10,  20 
and  40  nanoseconds  (ns)  and  by  comparing  the  FFT  of  the  resulting  data. 
The  results  are  shown  in  Fig.  10.  The  response  at  10  and  20  ns  is  very 
close  but  at  40  ns  loss  in  amplitude  and  change  in  frequency  content  is 
observed.  This  shows  that  for  the  frequency  range  of  interest  in  this 
work  ('  10  MHz)  digitizing  intervals  of  10  or  20  ns  can  be  safely  used. 


fcN  k***  . 


EFFECT  OF  SAMPL/A/O  INTERVAL  ON 

TRANSDUCER  RESPONSE 


l-ii;.  10.  Effect  of  sampling  interval  on  transducer  response 


In  all  the  measurements  made  to  date  digitizing  intervals  of  10  or  20  ns 


have  been  used. 

The  second  factor  considered  in  analyzing  data  is  the  resolution  of 
the  signal  in  frequency  domain.  In  the  Fourier  transformation  the 
difference  in  frequency  between  the  digitizing  points  is  Zif=l/L  where  L 
is  the  signal  length.  The  signal  length  can  be  computed  by  L=NAt  where 
N  is  the  total  number  of  digitizing  points  and  At  is  the  digitizing 
interval.  As  described  earlier,  the  digitizing  interval  of  10  and  20  ns 
is  being  used  and  is  therefore,  fixed.  Hence  we  can  increase  L  by 
increasing  N  and  get  a  better  resolution  i.e.  smaller  Af  in  the 
frequency  domain.  A  resoltuion  of  Af^O.!  MHz  was  considered  satisfactory 
and  could  be  obtained  by  taking  N='102i1  and  At=10ns. 

The  third  factor  is  the  determination  of  the  useful  frequency  range 
of  transducer.  The  FFT  of  the  time  domain  signal  f(t)  in  Fig.  5.  is 
shown  in  Fig.  11.  This  is  the  plot  of  |F*(oo)|  .  Note  also  that  F*(w) 
appears  in  the  denominator  of  equation  (60)  and  hence  when  the  values  of 
I  F*(to) I  are  very  small,  will  give  rise  to  erroneous  results.  Thus,  as 
seen  in  Fig.  11,  the  useful  range  of  frequencies  over  which  the 
calculations  of  equation  (60)  can  be  performed,  is  limited.  We  chose 
the  criterion  to  limit  |F*(w)|  to  25$  of  its  peak  value.  The  range  so 
obtained  for  the  10MHz  transducer  on  an  aluminum  specimen  is  indicated 
in  Fig.  11.  It  is  to  be  noted  here  that  the  response  of  the  transducer 
has  some  meaning  when  it  is  being  used  against  a  specimen.  Hence  the 
same  transducer  will  respond  differently  against  different  specimens  and 
the  useful  range  has  to  be  obtained  for  every  transducer-specimen 
combination. 
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The  fourth  faotor  considered  is:  what  do  wo  moon  when  we  say  first 
or  second  pulse?  In  a  continuous  signal  where  do  we  decide  the  first 
pulse  has  ended  and  second  started?  The  answer  to  this  problem  is  nvt 
straightforward  and  especially  in  thin  samples  and  composite  materials 
there  is  no  logical  way  of  deciding  this.  Hence,  in  the  technique 
development  stage  we  took  aluminuiii  samples  of  thickness  such  that  one 
pulse  dies  down  completely  before  the  second  arrives.  Then  the  middle 
point  between  the  pulses,  where  the  signal  is  small,  is  taken  as  the  end 
of  one  pulse  and  the  start  of  the  next. 

This  technique  was  applied  to  aluminum  samples.  Tests  were 
conducted  and  wave  speeds  repeatable  to  +  0.05?  were  obtained.  One  very 
obvious  source  of  error  was  the  oil  layer.  The  thickness  of  the  oil 
layer  could  not  be  controlled.  This  source  of  error  was  eliminated  by 
immersing  the  specimen  and  transducer  in  water  and  using  water  as  the 
coupling  medium.  The  repeatability  of  the  tests  was  improved  from 
±  0.05?  to  ±  0.02?. 

In  the  next  phase  of  development,  the  technique  was  applied  to 
composite  specimen.  Following  problems  were  encountered 

a.  Composites  highly  attenuate  the  signals  being  used.  Hence 
the  second  pulse  (Fig.  12)  is  very  small  and  can  not  be  used  to  get  good 
results . 

b.  Composites  arc  highly  dispersive  and  hence  as  the  signal 
passes  through  the  specimen, it,  spreads  in  time  domain.  Due  to  this 
effect,  the  tail  of  one  pulse  and  the  head  of  the  subsequent  pulse  begin 
10  t  vc'rlap.  This  phenomenon  will  bo  cal''ed  the  spillover  of  the  pulses. 

The  first  problem  was  solv'^d  by  noting  from  equation  (50)  that  the 
front  .^urf  jce  r't  f  i  CL’t  i  on  fRay  2,  Fig.  2),  as  shown  in  Fig.  12.,  can  also 
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be  used  as  the  reference  wave  in  place  of  the  first  pulse.  Therefore, 
for  composites  we  use  the  front  surface  reflection  and  the  first  pulse 
to  calculate  the  wave  speed.  Since  front  surface  reflection  docs  not 
traverse  the  specimen,  it  does  not  suffer  from  dispersion  or  broadening. 
It  was  found  that  the  spillover  problem  was  minimized  when  the  front 
surface  reflection  and  the  first  back  surface  reflection  (first  pulse) 
were  used.  Hence  we  were  able  to  eliminate  problem  (a)  and  reduce 
problem  (b)  by  the  proper  selection  of  pulses. 

Two  major  problems  in  the  implimentation  of  this  are  being  resolved 
now;  these  are: 

(1)  The  operator  decides  when  one  pulse  ends  and  the  second  one 
starts.  This  slows  down  the  process  of  data  collection.  This  problem 
is  severe  and  is  unacceptable  for  thin  specimen. 

(2)  For  dispersive  medium  the  spillover  effect  becomes  very 
profound. 

Work  is  in  progress  aimed  at  eliminating  these  restrictions.  As  per 
equation  (50)  these  problems  can  be  eliminated  by  using  the  front 
surface  reflection  as  the  reference  signal  and  then  using  the  entire 
signal  with  all  the  reflections,  to  obtain  the  wavespeed.  We  are  also 
in  the  process  of  calculating  the  attenuation  of  the  waves  through  the 
specimen  which  is  another  important  factor  in  assessing  damage  in 


composites . 


.  Results* 


Figure  13  is  in  overview  of  the  equipment  used.  The  experiment  is 
initiated  at  the  pulser/recei ver  (Panametrics  Ultrasonic  Analyzer  Model 
5052UA)  The  analyzer  is  a  broadband  ultrasonic  device  which  includes  a 
pulser,  a  receiver,  a  stepless  gate,  and  a  gated  peak  detector.  The 
analyzer  sends  a  broadband  pulse  of  about  300  volts  peak-to-peak  to  the 
transducer  at  a  repetition  rate  selected  by  the  operator.  At  the  same 
instance  (taken  as  time  t=o)  a  triggering  pulse  is  sent  to  the 
digitizing  oscilloscope  (Data  Precision  Model  data  6000).  The 
transducer  launches  a  compressional  wave  in  water  towards  the  specimen. 
The  same  transucer  is  used  to  receive  the  reflected  signals  from  the 
specimen.  The  signals  are  then  sent  back  through  the  analyzer  where  the 
stepless  gate  is  used  to  isolate  the  desired  portion  of  the  waveform 
from  the  spurious  signals.  From  the  analyzer  the  gated  waveform  is  sent 
to  the  Data  6000  for  signal  processing.  All  FFT  computations  are 
performed  at  the  Data  6000  and  the  results  are  sent  to  the  Mine  PDP 
11/23  through  an  IEEE  488  bus  line. 

Longitudinal  wave  speed  in  the  composite  was  measured  by  using  the 
front-surface  reflection  (Ray  2,  Fig,  2)  and  the  first  back-surface 
reflection  (Ray  6,  Fig.  2)  and  eq.  (39).  The  initial  parameters  that 
must  be  determined  prior  to  running  the  program  are;  1)  frequency  range 
for  the  least  squares  curve  fit  of  the  phase  vs.  frequency  plot,  and  2) 
cursor  location  i.e.  the  point  at  which  the  two  reflections  will  be 
separated . 

All  composite  specimens  were  fabricated  using  Magnamite  AS4/35C.? 
prepreg  tape,  A  list  of  specimens  and  associated  geometry  is  given  in 
Table  1 . 

*Propared  by  .l.C..  Eden  (M.S.  Student) 


Tht  main  objective  for  using  FFT  signal  processing  was  to  achieve 
a  very  high  level  of  precision  in  wave  speed  measurement  (±0.2^).  To 
date  we  have  identified  four  major  sources  of  error:  1)  temperature,  2) 
water  absortion,  3)  positioning  of  specimen,  and  4)  repetition  rate. 

1.  Temperature  -  It  is  well  known  that  in  most  materials  the  stiffness, 
and  therefore  the  wavespeed  depends  upon  temperature.  In  order  to  carry 
out  a  well  "controlled"  experiment  it  was  deemed  essential  to  study  the 
temperature  dependence  of  the  wave  speed  in  a  composite  specimen.  The 
temperature  range  studied  was  from  6®C  to  30®C.  Specimen  16-5  was  used 
for  the  investigation  (see  Table  1).  From  Figure  14  it  is  observed  that 
over  the  range  of  6°  to  30°C  there  was  an  approximately  2%  decrease  in 
the  wavespeed.  The  straight  line  is  the  result  of  a  linear  least  square 
fit.  In  the  preseht  set-up  we  have  the  capability  to  control  the 
temperature  to  ±0.5°C.  From  Fig.  3  it  was  found  that  a  0.5®C  change  in 
temperature  corresponded  to  only  0.05?  change  in  velocity  i.e.  the 
temperature  contribution  to  the  total  error  is  about  0.05?  which  -  in 
view  of  the  overall  objective  of  ^  0.2?  error  -  is  an  acceptable  figure. 

2.  Water  Absorption  -  For  [0/90]  geometry  laminates  the  major  form  of 
damage  created  by  tensile  loading  is  transverse  cracking.  Since  this 
technique  involves  water  immersion  it  was  of  some  concern  that  the 
damaged  composite  may  be  absorbing  moisture  via  the  newly  created 
cracks.  It  was  found  that  for  highly  damaged  specimens  water  absorption 
during  measurements  could  cause  variations  of  _+3?  (of.  measurement  error 
of  0.2?).  In  order  to  eliminate  this  problem  specimen  edges  were  coated 
with  a  strippable  coating.  This  precaution  measurably  eliminated  the 
problem . 
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TEMPERATURE  INFLUENCE  TEST  \  16-5 


wavcsnccd  decreases  with  temperature 


Positioning 


Due  to  the  inconsist..-noie3  within  a  composite  (ie. 


3  • 

matrix  or  fiber  rich  regions,  thickness  variations,  surface  mark,  etc.) 
a  means  for  reproducing  the  same  position  of  the  transducer  relative  to 
the  specimen  was  developed.  It  is  estimated  that  the  position  accuracy 
iswithinj^O.OOlin. 

Repetition  Rate  -  As  mentioned  earlier  a  repetition  rate  is  selected 
by  the  operator.  In  the  early  stages  of  this  work,  we  selected  a 
repetition  rate  of  5KHz,  for  no  particular  reason.  This  means  that  the 
experiment  was  repeated  every  0.2  milliseconds  (ms).  Random  errors  as 
high  as  +5*  were  observed.  A  careful  examination  of  the  experimental 
procedure  revealed  the  following  source  of  the  error.  Suppose  the 
transducer  is  energized  as  time  t=0.  It  will  take  some  time  (say  t^) 
before  the  mechanical  energy  -  propagating  in  the  form  of  waves  -  will 
be  completely  dissipated.  If  the  subsequent  experiment  is  initiated  at 
some  time  t<tQ,  the  results  will  be  erroneous:  the  waves  leftover  from 
the  "previous"  experiment  will  interfere  with  the  waves  of  the  "present" 
experiment.  This  phenomenon  sets  an  upper  limit  on  the  repetition  rate 
(lower  limit  on  the  time  between  two  experiments).  By  trial  and  error  we 
found  the  upper  limit  on  repetition  rate  to  be  1  KHz  (i.e.  time  interval 
of  one  millisecond). 

Technique  Calibration  In  order  to  calibrate  the  technique  the 
longitudinal  wavespeed  of  aluminum  was  measured  first.  The  surfaces  of 
a  block  of  aluminum  (0.1895  in.  thick)  were  polished  to  a  0.1  pm  finish. 
The  water  bath  temperature  was  held  fixed  to  +_0.1°c.  A  transducer  with 
a  5MHz  center  frequency  was  used.  The  front-surface  and  the  first  back- 
surface  reflections  are  shown  in  Fig.  15;  the  FFT  of  the  former  is  shown 
in  Fig.  1o.  The  . specimen  war,  held  fixed  at  one  place.  The  results  of 


Fig.  15.  Front-surface  and  back-surface  reflections  from  an  aluminum  specimen 


five  measurementa  are  shown  in  Fig.  17;  the  reproducibility  is  ^0.02% 
i.e.  two  parts  per  ten  thousand.  Solid  lines  indicate  the  error  bounds. 
Next,  we  studied  the  effect  of  removing,  drying,  and  replacing  the 
specimen  (mimicing  the  corresponding  events  for  an  actual  composite 

specimen)  The  results  are  shown  in  Fig.  18.  The  precision  is  reduced 
from  0.02^  in  Fig.  17  (specimen  not  removed)  to  0.08?  here,  the  hash 
marks  indicate  the  error  bounds. 

The  attention  is  now  turned  to  a  Graphite/Epoxy  specimen  (16-5). 
To  study  the  "worst  case"  we  took  the  most  severely  damaged  specimen  . 
Fig.  19  shows  the  time  domain  signal  for  the  first  two  reflections. 

Note  that  due  to  dispersion  of  the  wave  the  second  pulse  has  "spilled 
over"  into  the  first  pulse.  The  operator  judgement  was  used  to  place 
the  cursor  mark  at  the  "end"  of  the  first  pulse.  The  FFT  of  the  first 

signal  is  shown  in  Fig.  20.  When  the  specimen  was  not  removed  from  its 

position,  the  repeatability  was  found  to  be  ±0.02?  (Fig.  21).  However, 
when  the  specimen  is  removed,  dried  and  replaced  the  error  increases  by 
a  factor  of  five  to  0.1?  (Fig.  22). 

In  conclusion,  the  present  technique  has  been  developed  to  the 
point  where  even  for  the  worst-damage  case  the  measurement  error  is 
claimed  to  be  about  0.2?  (using  a  factor  of  safety  of  two). 

Work  in  Progress 

Work  is  now  in  progress  in  which  a  specimen  will  be  loaded 
monotonically  to  a  certain  load  level  and  then  unloaded.  Following 
measurements  will  be  made:  (1)  Longitudinal  wavespeed; 

(2)  Axial  stiffness 

(3)  Edge  replication  (crack  density  and 
total  crack  surface  area) 


SECONDS 

reflections  from  a  highly  damaged  composite  specimen 


TEST  NUMBER 


The  specimen  is  then  reloaded  to  a  higher  load  level  and  the  sequence  is 


repeated  until  the  final  failure  occurs. 

Some  preliminary  results  have  been  obtained.  But  these  are 
contrary  to  what  was  expected.  These  results  are  being  carefully 
examined  for  some  additional  sources  of  error. 


5.  Professional  Personnel 


1.  Dr.  Vikram  K.  Kinra,  Principal  Investigator 

2.  Mr.  J.  Greg  Eden,  M.S.  Student 

3.  Mr.  Vinay  Dayal,  Ph.D.  Student 

4.  Mr.  Vasu  Iyer,  Ph.D.  Student  (Will  replace  Mr.  Eden  after  he 

finishes  his  M.S.  in  May/June  1985) 

5.  Mr.  J.  Grille,  Laboratory  Technician-Mechanical 

6.  Mr.  P.  Patti,  Laboratory  Technician  Electronics 
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1.  Attended  "Vibration  Damping  Workohop",  Long  Beach,  California 
February  1984. 

2.  Chaired  a  session,  "Composites  11"  at  the  V  International  Congress 
on  Experimental  Mechanics,  Montreal,  June  1984. 

3.  Presented  a  paper'  at  the  Sixth  International  Conference  on  Fracture 
Mechanics,  New  Delhi,  India,  December  1984. 

4.  Organiser  and  Chairman  of  a  session,  "Ultrasonic  Nondestructive 
Evaluation  of  Composites"  at  the  forthcoming  Spring  conference  of 
Society  for  Experimental  Mechanics,  Las  Vega,  June,  1985. 

5.  Will  be  presenting  an  invited  review  paper,  "Some  Experiments 
Concerning  Multiple  Scattering  of  Elastic  Waves  in  Random  Media"  at  a 
Symposium/Workshop,  "Multiple  Scattering  of  Waves  in  Random  Media  and  by 
Random  Rough  Surfaces",  The  Pennsylvania  State  University,  July  1985. 

6.  Held  technical  consultations  concerning  the  project  with  technical 
staff  of  the  following; 

(a)  Martin  Marietta  Aerospace,  Denver 

(b)  Douglas  Aircraft  Corporation,  Long  Beach 

(c)  Hercules  Corporation,  Utah 
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Nondestructive  Testing  of  Composite  Materials 


by 


John  Gregory  Eden 


The  methods  for  detecting  damage  may  be  separated  into 
two  categories:  1)  Destructive  testing,  and  2)  Nondestruc¬ 
tive  Testing  and  Evaluation  (NOTE).  Destructive  testing  is 
self-explanatory.  One  such  destructive  test  is  the  deply 
technique  [^].  Composite  test  specimens  are  pyrolized  after 
being  impregnated  with  gold  chloride  diethylether  solution. 
Upon  unstacking,  the  laminate  damage  due  to  loading  may  be 
viewed  with  the  use  of  a  stereo  microscope  and  flourescent 
lamp.  Destructive  testing  has  one  major  drawback-  the 
specimens  are  destroyed.  Therefore,  NDTE  techniques  are 
much  more  desirable  from  the  veiwpoint  of  field 
appl ications. 

Reifsnider,  et.al.  [5]  defines  NDE  as  '...that 
activity  associated  with  experimental  schemes  used  to 
interrogate  the  state  of  stress  and  state  of  the  material 
without  influencing  either  of  those  states  or  altering  the 
strength,  stiffness,  and  life  of  the  laminates  being 
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evaluated*.  A  complete  description  and  comparison  of  NDE 
techniques  is  given  by  Sendeckyj  [6].  Nondestrut  i  ve 
techniques  such  as  surface  replication,  X-ray  radiograghy, 
V i br o t her mo gr aghy ,  holograghy,  and  liquid  crystal  coatings 
have  been  applied  to  composite  materials  by  Stalnaker  [7], 
Dance  [8],  Jones  [9],  Maddux  [10],  and  Charles  [11], 
respectively. 

The  main  focus  of  the  present  researh  concerns  the 
application  of  ultrasonics  to  assess  damage  in  composite 
materials  due  to  fatigue  or  monotonic  loading.  Much  work 
has  been  done  in  the  area  of  Acoustic  Emission  (AE). 
Whenever  a  failure  mechanism  takes  place  (ie.,  fiber 
breakage,  matrix  cracking,  etc.)  stress  waves  are  emitted 
from  the  source.  By  monitoring  a  test  specimen  under 
loading  with  a  triangular  array  of  transducers  one  may 
locate  the  failure  site  by  using  the  arrival  times  of  the 
stress  pulses  received  at  the  transducers.  Bailey,  Freeman 
and  Hamilton  [12]  used  the  above  technique  in  conjuction 
with  X-ray  radiography  to  correlate  AE  signal  amplitudes 
with  different  types  of  damage  progression  in  fatigue 
loaded  composite  plates.  R.  Williams  [13]  has  combined  AE 
and  thermography  to  monitor  damage  growth  in  Boron/epoxy 
and  Bor  on/ al um i n urn  specimens  subjected  to  fatigue  loading. 
Arora  and  Tangri  [1^]  used  AE  techniques  to  detect  and 
continuously  monitor  subcritical  crack  growth  in  Zr-2.5? 
Nb.  Block  [15]  tested  a  number  of  composite  specimens  and 


determined  that  AE  events  generated  by  fiber  fracture  havj 


substantially  higher  peak  amplitudes  than  those  generated 
by  resin-controlled  mechanisms.  Ulman  [16]  observed  that  an 
increase  in  AE  count  rate  coincided  with  the  onset  of 
matrix  plasticity  and  also  damage  developement  at  high 
loads  in  metal  matrix  composites. 

As  reported  by  Williams  and  Egan  [17]  AE  results 
contained  in  much  of  the  literature  are  largely  qualitative 
even  though  quantitative  measures  are  often  pr^esented. 
Comparison  of  AE  research  is  almost  impossible  due  to  the 
insufficient  reporting  of  the  technique  used.  Also,  due  to 
the  complexity  of  the  problem,  even  within  a  single 
specimen  a  small  change  in  stress  may  result  in  a  number  of 
AE  events,  some  of  which  may  initiate  from  different 
sources.  Therefore,  comparison  of  single  AE  events  to  one 
another  will  most  likely  not  lead  to  source  mechanism 
discrimination.  However,  by  combining  individual  AE 
spectral  densities  to  derive  'mean'  normalized  densities 
which  in  turn  could  be  statistically  analyzed,  Williams  and 
Egan  were  able  to  provide  quantitative  discrimination 
between  AE  from  10,  90,  [±45,±45]  degree  specimens  during 
tensi le  loading . 

Another  type  of  ultrasonic  evaluation  is  the  C-scan. 
Blake  [18-21]  has  reported  on  the  application  of  C-scan 
digital  systems  to  NDE  of  composite  materials.  Analyzing 
the  signature  formed  by  digitally  monitoring  the  amplitude 


of  a  particular  segment  of  a  received  ultrasonic  waveform 
can  lead  to  information  regarding  material  discontinuities. 
Time  domain  representations  of  the  wave  form  can  be  used  to 
correlate  discontinuities  which  appear  similar  in  the 
amplitude  signature.  The  technique  involves  coupling  a 
transducer  to  the  specimen  and  propagating  longitudinal 
waves  through  the  thickness  direction.  The  same  transducer 
is  used  to  receive  the  reflected  sound  waves.  The  received 
signal  may  now  be  processed  by  analog  or  digiCally.  If  by 
analog,  the  gated  portion  of  the  total  waveform  is 
peak-height  analyzed.  The  peak  output  is  quantitized  to  ten 
discrete  levels  in  0.1  volt  increments.  Each  increasing 
level  corresponds  to  a  darker  shade  of  gray  when 
transmitted  to  a  pen  amplifier.  The  pen  burns  the  surface 
layers  of  an  i nk- i mpr egnat ed ,  electrically  conductive  paper 
to  produce  a  picture  of  the  specimen's  material  integrity. 
Alternatively,  one  could  digitize  the  received  signal 
[18-21].  The  digitized  signal  can  then  be  processed  by  a 
computer.  Peak  data  (voltage)  is  sent  through  a 
microprocessor  which  performs  gray  scale  conversions  and  in 
turn  sends  this  data  to  either  a  dot-matrix  printer,  an 
analog  conductive  paper  system,  or  a  color  graphics  system. 
Chang,  et.al.  [22,231  have  developed  an  in-service 
inspection  system  (ISIS)  which  can  produce  hard-copy  real 
time  plots  and  post-inspection  C-scans  to  be  used  in 
production  hand  scanning. Post-inspection  flaw 


magnification,  f 1  aw- ampl i t ude  listing,  and  RF  (radio 


f requency ) waveform  digitization  are  the  major  advantages  of 
the  ISIS.  Kiraly  and  Meyn  [2ii]  have  developed  a  computer 
controlled  scanning  system  to  monitor  the  initiation  and 
progression  of  local  damage  patterns  in  composite  specimens 
under  tensile  loading.  While  scanning  a  specimen,  a 
16-level  gray  scale  image  is  displayed  on  a  CRT  and  stored 
on  a  floppy  disk.  The  images  may  be  stored  in  sequence  of 
load  step  which  can  be  played  back  to  create  a  movie 

ft 

showing  damage  growth  as  a  function  of  loading.  Daniel  [25] 
has  conductea  C-scan  evaluations  of  fatigue  loaded 
specimens  with  temperature  and  moisture  variations  to 
monitor  flaw  growth  in  graphite/epoxy  laminates  containing 
initial  flaws.  Liber  [26]  has  used  C-scan  to  determine 
damage  behavior  in  flat  and  cylindrical  composite  specimens 
with  pre-existing  flaws.  The  limitation  of  the  C-scan 
technique  is  that  it  only  works  with  relatively  large 
delamination  flaws;  it  yields  little  information  about 
other  flaw  types.  However,  the  technique  is  qualified  in 
that  it  provides  a  pictorial  view  of  specimen  quality. 

In  the  work  done  by  Dreumal  and  Speijer  [27,28]  a 
technique  called  Polar  Scan  was  developed  which  can 
generate  an  image  that  is  a  unique  fingerprint  of  a 
specific  laminate.  The  technique  involves  rotating  the 
particular  specimen  of  interest  between  two  transducers  and 
monitoring  the  amplitude  of  the  through  transmitted 
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ringdown  oscillations  exceeding  a  preset  threshold  voltage 
Vary  has  shown  that  microvoid  content,  ultimate  tensile 
strength,  cure  pressure,  and  interlaminar  shear  strength 
all  effect  the  value  of  SWF.  Also,  the  SWF  was  determined 
at  various  positions  along  specimens  and  the  location  of 
minimum  values  of  SWF  correlated  with  the  actual  failure 
sites  when  the  specimens  were  loaded  to  failure. 

Following  Vary,  Williams  et.al.  [38-42]  have  studied 
the  effect  of  fatigue  and  impact  loading  on  SWF  and 
attenuation.  In  earlier  works  with  metals,  Truell  and 
Hikata  [43]  established  a  relation  between  the  attenuation 
and  number  of  recorded  fatigue  cycles  in  aluminum.  Joshi 
and  Green  [44]  have  successfully  used  attenuation  to 
monitor  fatigue  damage  in  pol y cr ys t al 1 i ne  aluminum  and 
steel  specimens.  Williams  and  Doll  [38]  have  monitored 
attenuation  and  wave  velocity  in  compression-compression 
fatigue  loaded  unidirectional  composites  in  the  frequency 
range  of  .5  to  2  MHz.  Within  ±5%  accuracy  no  change  in 
wave  velocity  as  a  function  of  fatigue  was  detected. 
However,  initial  attenuation  above  1.5  MHz  was  found  to  be 
a  good  indicator  of  relative  fatigue  life.  Williams  and 
Lampert  [42]  determined  that  impact  damage  in  graphite 
fiber  composites  can  be  assessed  quantitatively  using 
either  the  through -the-thickness  attenuation  or  the  SWF. 
Williams,  Hashemi,  and  Lee  [41]  conducted  a  complete  study 
of  ultrasonic  attenuation  and  velocity  in  AS/3501-6 


composites.  In  the  frequency  range  of  .25  to  MHz  they 


found  the  wave  velocity  to  be  frequency  independent 
although  attenuation  was  frequency  dependent.  Williams, 
Yuce,  and  Lee  [39]  measured  the  attenuation  at  ^  MHz 
through  composite  specimens  which  varied  in  cure 
temperature  and  pressures  and  found  good  correlation  with 
the  number  of  fatigue  cycles  to  failure  and  an  increase  in 
attenuation.  Hayford,  et.al.  [^5]  determined  correlations 
between  initial  attenuation  and  shear  strength  in  graphite 

ft 

polyimide  composites.  Henneke,  et.al.  [46]  have  also 
studied  the  relation  between  SWF  and  fatigue  life  in 
composite  materials.  Talreja  [47]  has  applied  spectral 
frequency  analysis  and  observed  close  agreement  between 
spectral  density  and  stiffness  changes.  Nimmer  [48]  has 
used  attenuation  and  wave  velocity  to  monitor  damage 
accumulation  in  composite  flywheel  disks.  Hemann  [49]  has 
conducted  a  study  of  the  effect  of  stress  on  ultrasonic 
pulses  in  fiberous  composites.  The  velocities  were  found 
to  be  frequency  dependent,  but  weakly  dependent  on  the 
stress  in  the  specimen. 

The  anisotropy  of  composite  materials  can  greatly 
complicate  the  interpretation  of  the  received  ultrasonic 
signal.  Kriz,  et.al.  [50-52]  and  Kinra  and  Eden  [53]  have 
applied  the  solutions  of  the  Christoffel  [54]  equations  to 
graphically  depict  stiffness,  longitudinal  wave  velocity, 
and  shear  wave  velocity  variations  in  specified  planes  for 
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graphite  epoxy  laminates.  Deviations  between  wave 
propagation  direction  and  energy  propagation  were  found  to 
be  as  large  as  60  degrees.  This  has  some  rather  interesting 
implications  concerning  the  interpretation  of  received 
ultrasonic  pulses.  Kriz  [51],  Tauchert  and  Guzelsu  [55,56] 
have  experimentally  determined  elastic  moduli  by  relating 
wave  velocity  measurements  in  various  propagation 
directions  to  stiffness. 

The  previous  velocity  measur emi^nts  in  composites  have 
relied  on  monitoring  the  time  delay  for  a  through 
transmitted  pulse.  This  technique  works  well  if  the 
specimen  is  thick.  If  too  thin,  then  multiple  reflections 
within  from  the  front  and  back  surfaces  may  interfere  with 
the  received  wave  packet.  In  order  to  alliviate  this 
problem  a  new  technique  ,  Spectroscopy,  has  been 
investigated.  Spectroscopy  utilizes  the  frequency  contents 
of  the  reflected  signals  instead  of  the  time  domain  data. 
Chang,  et.al.,  have  utilized  the  connection  of  resonance 
and  frequency  to  determine  wave  velocities  in  thin 
composite  specimens  within  ±2%  [57].  Since  Spectroscopy 
techniques  reley  on  signal  processing  of  digitized 
waveforms,  which  has  only  recently  been  established,  the 
amount  of  published  articles  is  few.  However,  much  work  is 
currently  being  done  in  this  area. 
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INTRODUCTION 


The  present  work  is  a  part  of  an  overall  effort  to  devise  experimental 
techniques  for  Nondestructive  Testing  and  Evaluation  of  fiber-reinforced  composite 
materials.  A  clear  understanding  of  the  nature  of  wave  propagation  in  anisotropic 
media  is  a  necessary  prerequisite  for  a  correct  interpretation  of  the  experimental 
data.  In  the  following,  we  explore  some  concepts  which  are  unique  to  composite 
materials:  velocity  surfaces,  slowness  surfaces,  and  deviations  of  group 

velocity  vector  and  particle  displacement  vector  from  the  wave  normal  vector. 


RESULTS  AND  DISCUSSIONS 


Tlie  anisotropic  nature  of  fiber-reinforced  composites  give  rise  to  some 
rather  interesting  behavior  patterns  when  plane  waves  are  propagated  througli 
them.  The  solutions  to  the  equations  of  Christoffel  [1]  are: 


where  C...  „  is  the  6x6  stifness  matrix,  n  a  unit  vector  normal  to  the  wavefront, 
ijki 

p  is  the  density,  v  is  the  phase  velocity,  6.,  is  the  Kronecker  delta,  and  p  is 

X  K 

a  unit  particle  displacement  vector.  These  equations  yield  three  characteristic 
eigen-velocities  with  which  are  associated  with  three  corresponding  eigea-vectors . 
As  the  wave  normal  takes  on  all  possible  values,  v(n)  vv/ill  trace  out  a  velocity 


surface  of  three  sheets.  A  computer  program  was  written  to  solve  the  above 
equations  for  the  most  general  composite  material.  The  numerical  values  presented 
here  are,  however,  only  for  the  case  of  a  unidirectional  composite  which  has  been 
modelled  as  a  transversely  isotropic  material.  The  five  independent  elastic 
constants  required  to  describe  such  a  material  are  presented  in  Table  1.  These 
are:  Young's  moduli!  and  E2»  shear  modulus  Poison's  ratios  v and 

v^^'  Figure  1  defines  all  angles  to  be  used  in  the  following  sections.  Figures 
2-6  show  plots  of  the  three  velocity  sheets  for  each  material  as  well  as  the 
particle  displacement  vector,  p.  Because  of  symmetrices  of  a  transversely 
isotropic  material,  only  a  quarter  of  the  intersection  of  the  velocity  surfaces 
with  the  plane  needs  to  be  plotted.  From  Figs.  2-6  it  is  obvious  that 

along  one  branch  the  phase  velocity  increases  greatly  in  the  direction  of  the 
fibers  (X^  denotes  the  fiber  direction  while  X^  is  normal  to  the  fibers).  Since 
on  this  surface  the  particle  displacement  at  G  =  0°,  90°  is  purely  longitudinal 
this  velocity  surface  will  be  dent>ted  as  the  longitudinal  (L)  velocity  sheet. 

The  velocity  surface  which  has  the  particles  moving  transversly  in  the 


plane  will  be  called  the  transverse  (or  shear)  vertical  (SV)  and  the  last  surface 


with  thu  particles  moving  transversly  in  the  X^-direction  will  be  referred  to  as 
the  transverse  (or  shear)  horizontal  (SH) .  As  expected  Figs.  2-6  indicate 
that  as  stiffness  in  the  1,  2,  or  3  direction  increases  or  decreases  the 
corresponding  velocity  in  those  directions  will  follow  the  same  trend.  It 
may  also  be  observed  that  for  the  stiffer  materials  the  p-vector  rotates  much 
faster  towards  the  direction  of  highest  stiffness  when  0  is  increased. 

The  next  set  of  figures  (7-11)  show  the  deviation  (a)  of  the  particle 
displacement  vector  from  what  it  would  be  in  the  case  of  an  isotropic  material. 
For  the  longitudinal  branch  the  p-vector  is  parallel  to  the  wave-normal  whereas 
for  the  SV  and  SH  branches  the  p-vector  is  perpind icular  to  the  wave-normal. 

These  figures  show  only  two  curves  each,  but  actually  threei  curves  exist.  The 
L  and  SV  deviations  are  exactly  identical  while  the  SH  indicates  no  deviation 
from  the  isotropic  case.  It  is  often  common  to  define  a  particular  branch  on 
the  basis  of  the  direction  of  particle  displacement.  That  is,  if  the  angle 
between  the  p-vector  and  wave  normal  is  between  0®  and  45®  the  branch  should 
be  called  quasi-longitudinal  and  if  the  angle  is  between  45®  and  90®  it  would 
be  called  quasi-transverse  (0°  =  purely  longitudinal,  90°  =  purely  transverse). 
Figures  7-11  indicate  that  these  labels  may  be  misleading  due  to  the  tendency 
for  u  to  increase  beyond  45°  for  the  L  branch  and  then  decrease.  It  is  for  the 
above  reason  that  the  previously  defined  L,  SV,  and  SH  labels  were  based  upon 
0  =  0°,  90°  p-vector  directions.  As  the  material  stiffness  increase  so  does 
the  deviation  angle  which  again  indicates  the  preference  of  particle  motion  in 
the  direction  of  greater stif f ness  for  the  L  branch  and  in  the  direction  of  higher 
shear  strength  for  the  SV  branch. 

Slowness  is  defined  as  the  inverse  of  the  velocity.  Figures  12-16  show 
plots  of  the  slowness  surfaces  for  each  composite.  Also  on  the  plots  are  the 
directions  of  the  energy  propagation  vector  F.  given  by  [1] 


I  ^  I  ^ 


rr®" 


1  2  2 

F.  =  c  w  G  .  ,  p  .p,  s 

1  2  ijk;.  j  k  t 

I  2  ’ 

without  loss  of  goner. il  ity  the  factor  —  c  w“  w.i.s  set  equal  to  one  for  our  case, 
since  c  and  w  are  constants,  *^he  stiffness  matrix,  p  is  the  particle  dis¬ 

placement  vector,  and  s  is  the  slowaiess  vector  defined  as  s^  =  n^/v.  As  expected 
the  normal  to  the  slowness  surface  coincides  with  the  direction  of  the  energy 
vector  [!]•  The  plots  also  show  that  for  the  L  branch  the  energy  velocity 
direction  moves  rapidly  towards  the  fiber  direction  as  the  wave  normal  angle 
(G)  is  increased  from  0°  to  20°.  This  interesting  phenomenon  is  shown  explicitly 
in  Figs.  17-21  where  the  deviation  angle  (A)  between  the  energy  vector  and  wave 
normal  is  plotted.  The  plots  also  show  an  interesting  feature  along  the  SV  branch. 
The  energy  initially  shows  a  preference  to  propagate  in  the  'direction  of  the 
fibers  for  small  values  of  0.  However,  as  0  is  further  increased  there  is  a  cross¬ 
over  point  beyond  which  the  energy  prefers  to  propagate  away  from  the  fiber-direction. 
At  this  time  we  have  no  reasonable  explanation  for  this  phenomenon. 


CONCLUSIONS 


In  conclusion,  the  dynamic  response  of  unidirectional  fiber-reinforced 
composites  has  been  studied  in  some  detail.  This  study  has  revealed  several 
interesting  features  which  are  peculiar  to  anisotropic  media.  For  example,  the 
deviation  between  the  energy  vector  and  the  wave  propagation  direction  may  be 
as  large  as  60°  in  some  cases.  This  has  some  obviously  rather  interesting 
implications  concerning  excitation  and  detection  of  elastic  disturbances  in 
composite  materials  using  piezoelectric  crystals.  This  study  is  a  prelude  to 
Ultrasonic  Nondestructive  Testing  of  composite  materials. 
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